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Abstract 


The  problem  addressed  in  this  memorandum  is  nonstationary  interference  suppression 
in  noise  radar  systems.  Towards  this  aim,  linear  time- frequency  (TF)  transforms, 
short  time  Fourier  transform  (STFT)  and  local  polynomial  Fourier  transform  (LPFT) 
are  used  as  a  means  of  signal  representation.  The  noise  radar  return  signal  is  a 
wideband  random  signal  occupying  the  whole  TF  plane,  while  the  interference  signal 
is  well  concentrated  in  the  TF  plane.  This  implies  that  the  filtering  of  the  received 
signal  can  be  performed  by  using  a  binary  mask  to  excise  only  a  portion  of  the  TF 
plane  corrupted  by  the  interference.  Simulations  carried  out  on  the  radar  return  signal 
corrupted  by  an  extremely  strong  nonstationary  interferences,  covering  the  same  time 
and  frequency  ranges.  Results  confirm  the  effectiveness  of  the  proposed  method. 


Resume 


Le  present  document  aborde  le  probleme  de  la  suppression  du  brouillage  non  stat.ion- 
naire  dans  les  systemes  radar  a  onde  de  bruit.  Dans  ce  contexte,  des  transformees 
lineaires  temps-frequence  (TF),  des  transformees  de  Fourier  fenetrees  (STFT)  et  des 
transformees  de  Fourier  polynomiales  locales  (LPFT)  sont  utilisees  pour  representer 
les  signaux.  L’echo  du  radar  a  onde  de  bruit  est  constitue  d’un  signal  aleatoire  a  large 
bande  occupant  tout  le  plan  TF,  alors  que  le  signal  de  brouillage  est  bien  concentre 
dans  le  plan  TF.  Le  filtrage  du  signal  regu  peut  ainsi  etre  effectue  a  l’aide  d’un  masque 
binaire  qui  extrait  seulement  une  partie  du  plan  TF  corrompu  par  le  brouillage.  Des 
simulations  ont  ete  effectuees  au  moyen  d’un  echo  radar  corrompu  par  du  brouillage 
non  stationnaire  extremement  intense,  couvrant  les  memes  plages  de  temps  et  de 
frequence.  Les  resultats  ont  confirme  l’efficacite  de  la  met  hock  ;  proposee. 
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Executive  summary 

Nonstationary  Interference  Excision  for  Noise  Radar 
Systems  using  Time-Frequency  based  Methods 

T.  Thayaparan,  M.  Dakovic,  L.  Stankovic;  DRDC  Ottawa  TM  2007-334;  Defence 
R&D  Canada  -  Ottawa;  December  2007. 

Background:  Effective  battlefield  radar  surveillance  and  protection  of  the  nation’s 
borders,  harbours,  air  space,  and  marine  assets  require  covert  and  interference-free 
operation.  This  includes  not  only  that  friendly  radar  signals  be  undetected  by  the 
enemy,  but  also  that  such  systems  be  immune  from  jamming  and  external  electro¬ 
magnetic  interference  (EMI).  While  many  low  probability  of  intercept  (LPI)  and  low 
probability  of  detection  (LPD)  waveforms  have  been  developed  and  tested  over  the 
past  years,  these  primarily  rely  upon  the  use  of  pseudorandom  transmit  waveforms. 
However,  with  the  increasing  storage  and  computational  capabilities  of  advanced  di¬ 
gital  signal  processors,  it  has  become  easier  for  the  intelligent  adversary  not  only  to 
detect,  characterize,  and  recognize  such  signals,  but  also  to  use  the  information  to 
jam  and  confuse  friendly  radar  systems.  Random  noise  radar  is  an  attractive  and 
viable  option  for  use  in  these  applications. 

Random  noise  radar  refers  to  techniques  and  applications  that  use  incoherent  noise 
as  the  probing  transmit  waveform.  One  of  the  major  advantages  of  using  noise  as  the 
transmit  signal  is  its  inherent  immunity  from  detection,  unintended  interference,  and 
hostile  jamming. 

Results:  The  problem  addressed  in  this  memorandum  is  nonstationary  interference 
suppression  in  noise  radar  systems.  Towards  this  aim,  linear  time- frequency  (TF) 
transforms,  short  time  Fourier  transform  (STFT)  and  local  polynomial  Fourier  trans¬ 
form  (LPFT)  are  used  as  a  means  of  signal  representation.  The  noise  radar  return 
signal  is  a  wideband  random  signal  occupying  the  whole  TF  plane,  while  the  inter¬ 
ference  signal  is  well  concentrated  in  the  TF  plane.  This  implies  that  the  filtering 
of  the  received  signal  can  be  performed  by  using  a  binary  mask  to  excise  only  a 
portion  of  the  TF  plane  corrupted  by  the  interference.  Simulations  carried  out  on 
the  radar  return  signal  corrupted  by  an  extremely  strong  nonstationary  interference 
sources,  covering  the  same  time  and  frequency  ranges.  Results  confirm  the  utility  of 
the  proposed  method. 

Significance:  Results  clearly  demonstrate  that  the  time-frequency  based  interfer¬ 
ence  suppression  method  can  significantly  improve  the  detection  performance  of  the 
noise  radar  systems  in  extremely  strong  nonstationary  interference  environments. 
This  interference  suppression  method  can  potentially  be  extended  to  other  relevant 
applications  such  as  over-the-horizon  radar  systems  (OTHR),  through- wall  imaging 
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systems,  phased  array  radar  systems,  global  positioning  satellite  (GPS)  receivers, 
global  navigation  satellite  systems  (GNSSs),  etc.  Another  important  application  of 
this  method  is  jammer  rejection  in  the  spread  spectrum  (SS)  communication  systems. 
Different  methods  have  been  proposed  for  rejection  or  mitigation  of  interferences  of 
this  kind,  in  order  to  improve  interference  immunity  of  SS  systems  and  provide  more 
reliable  receiving  and  decoding  of  the  useful  signal.  The  proposed  method  presents 
another  viable  approach  to  jammer  mitigation  and  enhances  the  performance  of  the 
SS  receiver  in  such  severe  interfering  environment.  The  proposed  method  may  also 
be  extended  to  the  case  of  multiple  jammers.  This  method  can  generally  be  applied 
whenever  the  desired  signal  is  corrupted  by  broadband  interferences  characterized  by 
narrowband  instantaneous  bandwidths. 


IV 
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Nonstationary  Interference  Excision  for  Noise  Radar 
Systems  using  Time-Frequency  based  Methods 

T.  Thayaparan,  M.  Dakovic,  L.  Stankovic;  DRDC  Ottawa  TM  2007-334;  R  &  D 

pour  la  defense  Canada  -  Ottawa;  decembre  2007. 

Contexte  :  Afin  de  surveiller  et  de  proteger  efficacement  les  frontieres,  les  ports, 
l’espace  aerien  et  les  ressources  maritimes  d’un  pays,  un  radar  de  champ  de  bataille 
doit  pouvoir  fonctionner  furtivement  et  sans  brouillage.  II  est  done  non  seulement 
necessaire  que  l’ennemi  ne  puisse  pas  detecter  les  signaux  du  radar  ami,  mais  aussi 
que  le  systeme  soit  protege  cont.re  Interference  et  le  brouillage  electromagnet ique 
externe.  Bien  qu’un  grand  nombre  de  formes  d’onde  a  faible  probability  d’interception 
(LPI)  et  faible  probability  de  detection  (LPD)  aient  ete  etablies  et  mises  a  l’essai  au 
cours  des  dernieres  annees,  elles  se  fondent  principalement  sur  des  formes  d’onde  a 
emission  pseudo-aleatoire.  L’accroissement  des  ressources  de  stockage  et  de  calcul 
des  processeurs  de  signaux  numeriques  perfectionnes  facilite  t.outefois  la  tache  de 
l’adversaire  intelligent  qui  veut  non  seulement  detecter,  caracteriser  et  reconnaitre 
les  signaux  transmis,  mais  aussi  utiliser  l’information  pour  brouiller  et  perturber 
les  systemes  radar  amis.  Le  radar  a  onde  de  bruit  aleatoire  constitue  une  solution 
attrayante  et  viable  pour  ces  applications. 

Le  radar  a  onde  de  bruit  aleatoire  fait  appel  a  des  techniques  et  applications  utilisant 
le  bruit  incoherent  comme  forme  d’onde  d’explorat.ion  et  de  transmission.  L’un  des 
principaux  avant.ages  de  ce  mode  de  fonctionnement  tient  a  l’immunite  inherente  a 
la  detection,  au  brouillage  involontaire  et  a  1’ interference  ennemie. 

Resultats  :  Le  present  document  aborde  le  probleme  de  la  suppression  du  brouillage 
non  stationnaire  dans  les  systemes  radar  a  onde  de  bruit.  Dans  ce  contexte,  des 
transformees  lineaires  temps-frequence  (TF),  des  transformees  de  Fourier  fenetrees 
(STFT)  et  des  transformees  de  Fourier  polynomiales  locales  (LPFT)  sont  utilisees 
pour  representer  les  signaux.  L’echo  du  radar  a  onde  de  bruit  est  constitue  d’un  signal 
aleatoire  a  large  bande  occupant  tout  le  plan  TF,  alors  que  le  signal  de  brouillage  est 
bien  concentre  dans  le  plan  TF.  Le  filtrage  du  signal  regu  peut  ainsi  etre  effectue  a 
l’aide  d’un  masque  binaire  qui  extrait  seulement  une  partie  du  plan  TF  corrompu  par 
le  brouillage.  Des  simulations  ont  ete  effectuees  au  moyen  d’un  echo  radar  corrompu 
par  des  sources  de  brouillage  non  stationnaire  extremement  intense,  couvrant  les 
memes  plages  de  temps  et  de  frequence.  Les  resultats  ont  continue  l’efficacity  de  la 
methode  proposee. 

Portee  :  Les  resultats  demontrent  clairement  que  la  methode  temps-frequence  de 
suppression  du  brouillage  peut  ameliorer  considerablement  le  rendement  de  detec- 
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tion  des  systemes  radar  a  onde  de  bruit  dans  des  conditions  de  brouillage  non  sta- 
tionnaire  extremement  intense.  Cette  methode  de  suppression  du  brouillage  pourrait 
eventuellement  etre  adaptee  a  d’autres  applications  pertinentes,  comme  les  systemes 
radar  transhorizon  (OTHR),  au  moyen  de  systemes  d’imagerie  a  travers  les  murs, 
de  systemes  radar  a  commande  de  phase,  de  recepteurs  du  systeme  de  positionne- 
ment  a  couverture  mondiale  (GPS),  de  systemes  mondiaux  de  navigation  par  satellite 
(GNSS),  etc.  Cette  methode  trouve  une  autre  application  importante  dans  le  rejet 
du  brouillage  a  l’interieur  des  systemes  de  communications  a  spectre  etale  (SS).  Dif- 
ferentes  methodes  ont  ete  proposees  pour  le  rejet  ou  l’attenuation  du  brouillage  de  ce 
type,  de  maniere  a  ameliorer  Fimmunite  au  brouillage  des  systemes  SS  et  a  permettre 
la  reception  et  le  decodage  plus  fiables  du  signal  utile.  La  methode  proposee  presente 
une  autre  solution  viable  a  l’attenuation  du  brouillage  et  ameliore  le  rendement  du 
recepteur  SS  dans  des  conditions  de  brouillage  intense.  Elle  peut  aussi  s’adapter  au 
cas  de  brouillage  multiple.  Cette  methode  peut  generalement  s’appliquer  chaque  fois 
que  le  signal  utile  est  corrompu  par  du  brouillage  a  large  bande  constit-ue  de  bandes 
instantanees  etroites. 


VI 
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1  Introduction 


The  term  "random  noise"  as  applied  to  radar  refers  to  techniques  and  applications 
that  use  incoherent  noise  as  the  probing  transmit  waveform.  Because  of  the  truly 
random  transmitting  signal,  noise  radars  have  many  advantages  over  conventional 
radars,  including  unambiguous  estimation  of  both  range  and  velocity,  high  immunity 
to  noise,  low  probability  of  intercept  (LPI),  high  electromagnetic  compatibility,  good 
electronic  counter  countermeasure  (ECCM)  capability,  good  counter  electronic  sup¬ 
port  measure  (CESM)  capability,  and  ideal  thumbtack  ambiguity  function  [1]- [12] . 

Over  the  past  few  years,  the  research  has  been  devoted  to  the  development  and 
implementation  of  random  noise  radar  by  various  research  groups  [4],  [7],  [8],  [9]. 
Recent  research  has  investigated  the  potential  use  of  noise  radar  for  ultrawideband 
SAR/ISAR  imaging,  Doppler  and  polarimetric  measurements,  collision  warning,  de¬ 
tection  of  buried  objects,  and  targets  obscured  by  foliage  [2],  [5],  [8]- [12].  Wide  band¬ 
width  provides  high  range  resolution,  and  an  extended  pulse  length  reduces  peak 
power.  The  non-periodic  waveform  suppresses  the  range  ambiguity  while  reducing 
both  the  probability  of  intercept  and  interference. 

Mutual  interference  and  low  probability  of  interception  capabilities  of  noise  radar 
were  evaluated  in  previous  studies.  The  results  show  that  noise  radars  are  unlikely  to 
interfere  with  other  noise  radar  systems  or  other  radar  systems  in  the  same  band.  It  is 
also  shown  that  in  a  variety  of  noisy  environments,  the  noise  radar  has  a  much  lower 
LPI  than  the  conventional  LFM  radar.  The  noise  radar’s  exceptional  performance 
in  the  above  evaluations  indicates  that  it  is  a  suitable  radar  system  for  a  variety  of 
applications  frequently  improving  upon  the  performance  of  conventional  systems  [13], 

[14]- 

In  this  memorandum  we  have  studied  the  influence  of  an  extremely  strong  determin¬ 
istic  broad-band  interference  (signal  to  interference  ratio  as  low  as  -40dB),  covering 
the  frequency  and  time  ranges  of  the  operating  noise  radar.  A  time- frequency  (TF) 
based  interference  suppression  technique  is  developed  and  is  based  on  the  property 
of  time-frequency  representations  to  localize  signals  in  the  TF  plane.  Two  TF  trans¬ 
forms,  the  short  time  Fourier  transform  (STFT)  and  the  local  polynomial  Fourier 
transform  (LPFT),  are  used.  More  precisely,  time- varying  filters  based  on  the  STFT 
and  LPFT  are  developed.  Since  the  random  noise  radar  signal  occupies  the  whole 
TF  plane,  while  the  interference  signal  is  a  broadband  signal  characterized  with  a 
narrow  instantaneous  bandwidth,  the  time-varying  filtering  is  performed,  via  binary 
mask,  which  removes  the  interference’s  TF  signature  without  significant  degrada¬ 
tion  of  the  radar  return  signal.  Moreover,  the  LPFT  based  receiver  outperforms  the 
STFT  based  receiver  since  it  optimally  concentrates  an  interference  source  in  the 
TF  plane.  In  numerical  illustrations  we  have  considered  two  types  of  interferences:  a 
broad-band  sinusoidally  modulated  signal  and  a  linearly  frequency  modulated  (LFM) 
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signal,  corresponding  to  a  strong  interfering  LFM  radar,  which  covers  the  same  time 
and  frequency  ranges  as  the  operating  noise  radar. 

The  theoretical  background,  including  the  STFT,  the  LPFT,  time  varying  filtering, 
and  correlation-based  noise  radar  principles,  is  given  in  Section  2.  Section  3  introduces 
two  methods  for  binary  mask  implementation.  The  proposed  TF  filtering  methods’ 
performances  are  evaluated  by  means  of  numerical  examples  in  Section  4.  It  has  been 
shown  that  the  noise  radar  performs  in  a  satisfactory  way,  even  with  very  strong 
nonstationary  broadband  interference. 


2 
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2  Theoretical  background 


In  this  section,  a  short  introduction  to  linear  TF  methods,  i.e. ,  STFT  and  LPFT, 
is  given.  Furthermore,  the  section  also  provides  a  brief  description  of  time-varying 
filtering  and  correlation  based  noise  radar  principles. 

The  baseband  received  signal  r(n)  comprises  three  sequences  as  follows: 

r  (n)  =  x  (n)  +  I  ( n )  +  £  (n)  (1) 

where  x(n)  is  a  noise  radar  signal  sequence  (complex  white  Gaussian  noise  sequence 
with  zero  mean  and  variance  a2. ) ,  l  (n)  is  an  interference  signal  sequence  and  £(n) 
is  a  complex  additive  white  Gaussian  noise  (AWGN)  sequence,  with  zero  mean  and 
variance  cr|  uncorrelated  with  x(n).  The  interference  is  assumed  to  be  a  nonstationary 
signal  characterized  by  a  narrowband  instantaneous  bandwidth  and  by  the  following 
expression: 

I  (n)  =  Aiej<Pl{n) 

where  yq  (n)  is  the  phase  and  Aj  is  the  magnitude  of  the  interference. 

Signal-to- noise  ratio  (SNR)  and  signal-to-interference  ratio  (SIR)  are  defined  in  the 
following  way: 


SNR  = 

20  log10  — 

(2) 

SIR  = 

20  logio 

(3) 

The  influence  of  the  interference  signal  on  the  desired  radar  signal  can  be  mitigated 
by  using  TF  methods.  Herein,  we  are  interested  only  in  linear  TF  methods,  that  allow 
the  perfect  reconstruction  (synthesis)  of  the  observed  signal.  Two  such  methods,  i.e., 
short-time  Fourier  transform  and  local  polynomial  Fourier  transform,  will  be  used  as 
a  means  of  interference  suppression  in  this  memorandum. 

2.1  Short-time  Fourier  transform 

The  STFT  of  the  signal  r  (n)  [16],  denoted  as  STFTr  (n,  k).  is  obtained  by  sliding 
the  window  function  w{m )  over  the  signal  r  (n)  and  implementing  the  DFT  on  the 
product  of  r  (n)  and  window  at  the  current  position,  i.e., 

N/  2-1 

STFTr  (n,k)  =  ^  r  (n  +  m)iv  (m)  e~jwmk 

m=—N/  2 

=  DFT  [r  (n  +  m)  iv  (to)]  (4) 
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where  N  is  the  number  of  frequency  bins  adopted  in  the  DFT  calculation.  The 
window  function  is  usually  real  and  symmetric  with  the  property  that  w(0)  =  1.  A 
simple  manipulation  of  (4)  gives  the  STFT  synthesis  equation 


r  (n) 


1 

N 


N- 1 


E 


STFTr  (n,  k ) . 


(5) 


This  relation  states  that  the  signal  r(n)  can  be  obtained  by  summing  its  STFT  values 
over  the  frequency  variable  for  the  fixed  instant  n. 


2.2  Local  polynomial  Fourier  transform 


The  LPFT  has  been  recently  introduced  in  the  TF  analysis  [17],  [18],  with  the  Mth 
order  discrete  form  of  the  LPFT  of  the  sequence  r(n)  being  defined  as 


N/  2-1 

LPFT™  (n,  k )  =  ^  r  (n+m)w(m)  e 

m=—N/2 


ni+l 


M 

3  Z  (i+i)!  _  j  27T  u 
i= 1  C  J  N 


i+ 1 


=  DFT  r(n+m)w(m)e  i=1 


1  E  ™+i) 


(6) 


where  w(m )  and  N  are  the  same  as  in  the  STFT  definition  and  ujt  is  the  ?'th  transform 
parameter.  The  relation  (6)  indicates  that  the  LPFT  of  the  received  signal  can  be 
calculated  analogously  to  the  STFT,  i.e.,  by  sliding  the  analysis  window  w(m)  over 
the  modulated  received  signal 


M 

-3  E 


rin+mte 


’’  (*+!)■ 


and  implementing  the  DFT  on  the  product  of  the  modulated  signal  and  window  at 
the  current  position. 


The  LPFT  parameters  ay  for  i  =  1, 2, ...,  M  are  calculated  so  as  to  optimally  concen¬ 
trate  the  signal  (i.e.,  interference  in  this  case)  in  the  TF  plane  for  a  given  analysis 
window.  Towards  this  goal,  an  order  adaptive  algorithm  is  developed  in  [17]  and  it 
is  shown  to  keep  calculation  complexity  at  a  relatively  low  level.  Furthermore,  it  is 
shown  that  the  second-order  LPFT  produces  results  almost  independent  of  the  para¬ 
meters  of  FM  interference  sources,  thus  preventing  the  need  for  a  time-consuming 
calculation  of  the  higher-order  LPFT. 


2.3  Time-varying  filtering:  Binary  mask 

The  spectrum  of  the  noise  radar  signal  is  flat,  while  the  interference  signal  occupies 
a  narrow  frequency  band  at  each  time  instant.  The  time-varying  filtering  described 
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in  [15]  can  be  easily  implemented  here.  The  interference  excision  is  performed  in  the 
TF  plane  by  removing  its  TF  signature  through  a  time-varying  filter.  This  filter  can 
be  implemented  as  a  binary  mask,  denoted  as  B,  which  is  a  function  defined  in  the 
following  way: 


B  (n,  k ) 


0,  interference  exists  in  (n,  k ) 
1,  otherwise. 


(7) 


Practically,  B  (n,  k )  will  equal  1  only  for  points  (n,  k )  of  the  TF  plane  where  an 
interference  power  can  be  neglected. 


The  synthesis  is  performed  on  the  masked  transform  to  recover  the  "jammer- free" 
received  signal  r'  (n)  as  follows: 

N- 1 

r'  (n)  =  —  STFTr  (■ n ,  k)  B  (n,  A:)  (8) 

k=0 

or 

AT-l 

r'  (n)  =  -y  LPFTtm  (n, k)B  ( n,  A) .  (9) 

fc=0 


2.4  Correlation  Receiver 

The  correlation  receiver  uses  the  principle  that  when  the  reference  signal,  delayed 
by  Tref,  is  correlated  with  the  actual  target  echo,  the  peak  value  of  the  correlation 
function  indicates  the  distance  to  the  target  (the  amount  of  time  delay  of  the  reference 
signal  is  also  a  measure  of  distance  to  the  target),  while  Doppler  filters,  following 
the  correlator, 1  output  target  velocity  [7].  In  this  method,  the  return  signal  from 
the  target  is  cross- correlated  with  a  time-delayed  replica  of  the  transmit  waveform. 
When  Tref  is  varied  a  strong  correlation  peak  is  obtained  for  Tref  =  T{).  which  gives 
an  estimate  of  the  target  range  r0  =  cT0/ 2. 

Let  us  consider  a  radar  emitting  a  time-limited  signal  x(t).  Denote  the  received  signal 
by  y(t).  Furthermore,  we  assume  that  a  single  point  scatterer  is  located  at  the  range 
ro  along  the  radar  line-of-sight  (LOS).  From  this  assumption,  the  received  signal  can 
be  written  as: 

y(t)  =  Aax(t  -  T0)  +  e(t)  (10) 

where  To  =  2 r0/c  is  the  round-trip  delay  caused  by  the  finite  speed  of  the  electro¬ 
magnetic  waves,  s(t)  is  an  undesired  part  of  the  received  signal  (noise  caused  by  the 
reflection  from  other  objects  along  the  LOS  and  possible  jamming  signals)  with  A„ 
denoting  target  reflectivity.  Without  loss  of  generality  we  will  assume  that  Aa  =  1. 
The  correlation  of  the  emitted  and  received  signal  can  be  written  as: 

rTint 

R(t)  =  /  y(t)x*(t  —  r)dt.  (11) 

Jo 
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where  Tint  is  the  integration  time.  In  the  noiseless  case,  the  maximum  value  of  |i?(r)| 
occurs  at  the  point  r  =  T0. 


Let  us  now  assume  that  x(t)  is  a  white  stationary  Gaussian  random  process  with 
autocorrelation  function  Rx x(r).  The  output  of  the  correlation  receiver  given  by  (11) 
is  also  a  random  process.  Let  us  analyze  the  expected  value  of  (11)  as: 


E\mr)\ 


r Tint 

e[  /  y(t)x*(t 

Jo 

r Tint 

/  E[y(t)x*(t 

Jo 


r)dt] 

r)]dt 


E[x(t  —  T0)x*(t  —  t)]  +  E[e(t)x*(t 

rTint 

Rxx(r  —  T0)dt  +  /  E[e(t)x*(t- 


-  r)]dt 
r)]dt 


(12) 


Jo  Jo 

If  the  emitted  signal  x(t)  and  the  noise  z{t)  are  independent  processes  then  the  second 
term  in  (12)  is  equal  to  zero  and  we  get: 


E[R(t)]  =  TintRxx(r  -  T0). 


(13) 


Since  the  autocorrelation  function’s  maximum  occurs  at  u  =  0  (R(t)  <  R( 0)),  the 
delay  T0  can  be  estimated  as  the  position  of  the  maximum.  Thus: 

T0  =  max\E[R(T)]\  (14) 

T 


Special  cases: 

•  Let  x(t)  be  the  white  stationary  Gaussian  random  process.  The  autocorrela¬ 
tion  function  is  Rxx(r )  =  Io5(t  —  r).  This  is  an  ideal  shape  since  E[R(t)}  = 
TintI0S(t  —  r),  and  its  maxima  are  well  defined  (only  one  point  is  different  from 
zero).  Note  that  signals  of  this  form  are  not  bandlimited  and  they  can  not  be 
used  in  practical  applications. 

•  Let  x{t)  be  the  bandlimited  white  stationary  Gaussian  random  process  with 
power  spectral  density  (PSD)  Sxx(f )  =  So  for  /0  —  B /2  <  f  <  f0  +  B/2  and 
Sxx  ( / )  =  0  otherwise.  The  autocorrelation  function  is  of  the  form: 

JWr)  =  (15) 

7 TT 

with  well  a  defined  maximum  at  r  =  0,  and  with  a  first  side  lobe  that  is  -Be¬ 
times  lower  than  the  main  lobe. 
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3  Binary  mask  implementation 


This  section  offers  two  methods  of  binary  filtering  mask  implementation  and,  through 
numerical  examples,  assess  their  performances,  both  in  the  STFT  and  LPFT  case.  In 
all  examples  presented  in  this  section,  the  length  of  the  received  sequence  is  L  =  2048, 
N  =  256  and  SNR  =  0  dB.  The  interference  is  assumed  to  be  a  sinusoidal  FM  signal 
characterized  by  SIR  =  —20  dB.  In  the  LPFT  calculation,  the  perfect  knowledge 
of  the  LPFT  parameters  is  assumed.  Furthermore,  only  the  first  and  second  order 
LPFT  will  be  herein  considered. 


In  order  to  assess  performances  of  the  proposed  filtering  methods,  i.e.,  to  estimate 
a  remaining  interference  power  compared  to  a  remaining  radar  signal  power  after  a 
binary  mask  implementation,  the  following  ratios  are  introduced: 


L  N 

E  E\STFTx(n,k)Bs(n,k)\2 

SIRs  =  101og10=^ - 

^  j:\STFTj  (n,k)B 3  (n,k)\2 

n= 1 k=l 

£  K 

E  El LPFT}(n,k)BL1  (n, k)\ 

sirli  -  ioiog10^n=i -  (is) 

E  E  \LPFT}  (n,  k)  BL1  (n,  fc)|2 

n= 1 k= 1 

£  K 

E  J2\LPFT2(n,k)BL2(n,k)f 

sirl2  =  ioioglo==rti - 

E  E  \LPFTf  (n,  k)  Bl2  (n,  fc)|2 


where  Bs(n,k),  BLi(n,k )  and  BL2(n,k )  respectively  represent  binary  masks  ob¬ 
tained  in  the  STFT,  the  first  and  second  order  LPFT  based  filtering  procedures. 


3.1  Type  I  binary  mask 

The  first  adopted  binary  mask  is  trivial,  i.e.,  it  is  assumed  to  excise  all  frequency 
bins  of  the  transform,  whether  corrupted  by  interference  or  not,  that  exceed  some 
threshold  value.  The  following  threshold  value  will  be  assumed  [19]: 

T\  =  E  [| STFTxH  {n,  k)\2]  +  2\Jvar  [| STFTx+i  (n,  k)\2}  (17) 

where  E  [■]  and  Var  [•]  respectively  represent  the  expectation  and  variance  operator. 
Clearly,  STFTX+ %  (n,  k)  represents  the  STFT  of  the  sum  x  (n)  +  ^  (n). 
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Frequency  Frequency 


Figure  1:  Sinusoidal  FM  interference  case.  First  row:  Type  I  binary  masks  for  (a) 
STFT,  (b)  LPFT 1  and  (c)  LPFT2  based  interference  excision.  Second  row:  Type  II 
binary  masks  for  (d)  STFT,  (e)  LPFT 1  and  (f)  LPFT2  based  interference  excision. 
Zero  values  are  shown  in  black. 
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Figure  2:  LFM  interference  case.  First  row:  Type  I  binary  masks  for  (a)  STFT  and 
(b)  LPFT1  based  interference  excision.  Second  row:  Type  If  binary  masks  for  (c) 
STFT  and  (d)  LPFT 1  based  interference  excision.  Zero  values  are  shown  in  black. 
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The  advantage  of  this  type  of  binary  mask  is  a  simple  hardware  realization.  However, 
its  drawback  is  the  removal  of  a  certain  number  of  frequency  bins  that  are  not  cor¬ 
rupted  by  interference.  Moreover,  the  strongest  frequency  components  of  the  radar 
signal  are  eliminated  in  this  manner. 

The  first  and  third  rows  of  Table  1  give  values  of  SIR  ratios  (Equation  16),  averaged 
over  100  realizations.  Binary  masks  for  one  realization  of  the  STFT,  the  first  and 
second  order  LPFT  based  filtering  in  the  sinusoidal  FM  interference  case,  denoted 
as  B g,  B(  ,  and  B£2,  are  depicted  in  Figures  1(a),  (b)  and  (c),  respectively.  Cor¬ 
responding  binary  masks  for  the  LFM  interference  case  are  depicted  in  Figures  2(a) 
and  2(b)  with  superscript  I  denoting  the  first  type  of  binary  mask.  In  addition,  the 
percentage  of  the  excised  frequency  bins  for  all  implemented  transforms,  averaged 
over  100  runs,  is  shown  in  the  first  and  third  rows  of  Table  2.  Needless  to  say,  the 
second  order  LPFT  has  not  been  calculated  into  the  LFM  case  since  the  first  order 
LPFT  completely  focuses  LFM  interference  in  a  narrow  band. 

3.2  Type  II  binary  mask 

The  second  adopted  binary  is  more  sophisticated  than  the  first  one  as  it  removes  only 
corrupted  frequency  bins  around  the  spectral  peak  at  the  current  time  instant.  It  is 
obtained  by  means  of  the  following  steps. 

Step  1.  Set  the  binary  mask  B  (n.  k)  to  all  ones  and  set  n  =  1. 

Step  2.  If  n  >  L  exit;  otherwise  detect  the  maximum  of  |  ST FTr  (n,  k )  |  (or  |  LPFT (n,  k )  | ) 
at  the  current  time  instant  n.  Let  the  frequency  index  of  the  maximum  be  Aq 
and  set  k2  =  k±. 

Step  3.  As  long  as  | STFTr  (n,  k2) \  >  \STFTr  (n,  k2  -  1) |  or  \STFTr  (n,  k2) |2  >  T2, 
set  B  (n,  k2)  =  0  and  k2  =  k2  —  1. 

Step  4.  Set  k2  =  ki  +  1.  As  long  as  \STFTr  (n,k2)\  >  \STFTr  (n,  k2  +  1)|  or 
| STFTr  (n,  k2)\2  >  T2,  set  B  (n,  k2)  =  0  and  k2  =  k2  +  1. 

Step  5.  Set  n  =  n  +  1  and  go  to  step  2. 

Here,  T2  represents  the  threshold  value  defined  as 

T2  =  E  [| STFTx+i  (n,  A;)|2]  +  ^Var  [| STFTx+i  (n,  fc)|2].  (18) 

The  advantage  of  the  binary  mask  defined  in  this  way  is  the  sophisticated  interference 
removal.  We  begin  with  a  position  of  an  interference  spectral  peak,  at  the  needed 
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Table  1:  SIR  for  Type  I  and  Type  II  Binary  Masks 


Interference 

Binary  mask 

STFT 

LPFT1 

LPFT2 

Sin  FM 

Type  I 

7.49  dB 

11.66  dB 

13.76  dB 

Sin  FM 

Type  II 

16.08  dB 

20.41  dB 

25.3  dB 

Lin  FM 

Type  I 

12.62  dB 

16.58  dB 

— 

Lin  FM 

Type  II 

24.1  dB 

30.88  dB 

— 

Table  2:  Excised  Bins  Percentage  for  Type  I  and  Type  II  Binary  Masks 


Interference 

Binary  mask 

STFT 

LPFT1 

LPFT2 

Sin  FM 

Type  I 

33.49% 

12.03% 

8.4% 

Sin  FM 

Type  II 

33.32% 

10.18% 

6.03% 

Lin  FM 

Type  I 

14.1% 

6.81% 

— 

Lin  FM 

Type  II 

12.12% 

4.11% 

— 

time  instant  n,  and  set  B  (n.  k)  =  0  for  all  frequency  bins  k  of  the  transform  that  are 
corrupted  by  the  interference.  This  procedure  is  performed  as  long  as  the  absolute 
value  of  the  current  frequency  bin  is  greater  than  the  adjacent  one  or  its  squared 
absolute  value  is  greater  than  the  adopted  threshold  value  defined  by  (18).  Fur¬ 
thermore,  frequency  bins  of  the  transform  that  are  not  corrupted  by  an  interference 
remain  intact.  The  drawback  of  this  binary  mask  is  a  rather  complicated  hardware 
realization. 

The  second  and  fourth  rows  of  Table  1  present  obtained  values  of  ratios  (16),  averaged 
over  100  realizations.  Binary  masks  for  one  realization  of  the  STFT,  the  first  and 
second  order  LPFT  based  filtering  in  the  sinusoidal  FM  interference  case,  denoted  as 
B{/,  and  B((2,  are  depicted  in  Figures  1(d),  1(e)  and  1(f),  respectively.  Corres¬ 
ponding  binary  masks  for  the  LFM  interference  case  are  depicted  in  Figures  2(c)  and 
2(d).  Superscript  II  denotes  the  second  type  of  binary  mask.  Again,  the  percentage 
of  the  excised  frequency  bins  for  all  the  three  implemented  transforms,  averaged  over 
100  runs,  is  presented  in  the  second  and  fourth  row  of  Table  2. 

3.3  Discussion 

Despite  the  fact  that  the  removed  area  is  approximately  the  same  (see  Table  2),  the 
second  type  of  binary  mask  is  characterized  by  a  significant  improvement  in  SIR 
performance  compared  to  the  first  type,  as  shown  by  results  given  in  Table  1.  The 
reason  for  such  behavior  is  the  fact  that  the  first  type  of  binary  mask  eliminates 
the  strongest  frequency  components  of  the  radar  signal  along  with  the  interference 
components. 

The  STFT  based  filtering  is  outperformed  by  the  first  order  LPFT  based  filtering, 
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however,  the  second  order  LPFT  produces  the  best  results,  since  the  corresponding 
excised  area  is  the  smallest  compared  to  the  other  two.  Furthermore,  as  can  be  seen 
from  Figures  1(f)  and  2(d),  the  number  of  excised  frequency  bins  is  approximately 
the  same  for  each  time  instant,  meaning  that  increasing  the  LPFT  order  would  not 
provide  a  significant  SIR  improvement. 

Monocomponent  interference  has  been  assumed  in  this  analysis.  If  an  interference 
is  a  multicomponent  signal,  type  I  binary  mask  and  STFT  can  be  applied  without 
any  modifications,  while  the  procedure  for  the  type  II  binary  mask  implementation 
or  LPFT  needs  to  be  modified  in  order  to  remove  all  interference  components. 
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4  Simulations 


Consider  a  noise  radar  operating  at  carrier  frequency  /0  =  10  GHz  with  bandwidth 
B  =  204.8  MHz  and  pulse  duration  of  Tr  =  10  /is  (or  2048  samples).  The  received 
signal  is  sampled  at  the  Nyquist  rate  Ts  =  1/B.  The  radar  waveform  is  a  complex 
Gaussian  random  signal  with  i.i.d  (independent  and  identically  distributed)  real  and 
imaginary  parts.  The  transmitted  signal  is  reflected  from  the  single  point  scatterer 
target  located  at  distance  r0  =  1  km.  Let  us  also  assume  that  the  received  signal  is 
corrupted  by  a  Gaussian  complex  noise  £(f)  and  interference  signal  1(f).  The  received 
signal  is  of  the  form 

r(t)  =x(t)  +  I(t)+£(t)  (19) 

where  x(t)  =  Axxt(t  —  td)  represents  a  signal  reflected  from  the  target,  characterized 
by  attenuation  Ax  and  time  delay  t,i  =  — .  Note  that  r(n)  defined  by  (1),  represents 
the  discrete  version  of  (19). 

The  interference  is  assumed  to  be  a  frequency  modulated  signal  of  the  form 

I(t)  =  A^At) 

Two  types  of  interference  are  analyzed: 

1)  Interference  with  sinusoidal  instantaneous  frequency  (for  example  jamming  signal) 
in  the  form: 

/|i"(i)  =  wr^  =  “gBsill(8,ry)'  (20) 

2)  Interference  with  a  linear  instantaneous  frequency,  corresponding  to  a  LFM  radar 
operating  at  the  same  frequency  band.  The  linear  FM  interference  is  periodic  with 
period  and  the  equation  (21)  is  valid  within  the  fundamental  period  only. 

/f«)  =  ^  =  ~2 B)  for  -  |  <  t  <  )  (21) 

The  interference  suppression  is  performed  by  using  the  proposed  TF  based  filtering 
of  the  received  signal. 

Four  cases  are  considered 

1.  No  interference  suppression  is  performed. 

2.  Interference  suppression  is  performed  by  using  STFT. 

3.  Interference  suppression  is  performed  by  using  LPFT  of  the  first  order. 


DRDC  Ottawa  TM  2007-334 


13 


No  filtering  performed 


STFT  based  filtering 


range  range 


Figure  3:  Output  of  the  correlation  receiver  for  the  single  target  located  at  1000  m 
range  with  SNR  =  —20  dB  and  SIR  =  —20  dB  when  no  filtering  is  performed  (upper 
left),  STFT  based  filtering  is  performed  (upper  right)  first  order  LPFT  filtering  (lower 
left)  and  second  order  LPFT  (lower  right).  Type  II  binary  mask  is  used. 

4.  Interference  suppression  is  performed  by  using  LPFT  of  the  second  order. 

For  the  interference  suppression,  both  type  I  and  type  II  binary  mask  are  used,  as 
described  in  Section  3.  Figure  3  presents  radar  outputs  in  all  the  considered  cases. 
Note  that  the  target  is  not  detected  when  no  interference  suppression  is  performed, 
and  the  LPFT  based  filtering  outperforms  the  STFT  based  approach. 

The  simulation  was  performed  over  1000  realizations  of  the  received  signal.  The 
probability  of  false  target  detection1  versus  SIR  for  SNR  =  — 20dB  is  calculated  and 
results  are  presented  in  Tables  3  and  5  for  type  I  binary  mask  and  in  Tables  4  and  6 
for  type  II  binary  mask. 

1  False  target  detection  occures  when  detected  maxima  position  of  the  radar  output  does  not 
coincide  with  true  target  position. 
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SIR 

No  filtering 

STFT 

LPFT1 

LPFT2 

— 20dB 
— 25dB 
— 30dB 

27.7% 

46.9% 

15.3% 

11.0% 

75.7% 

57.1% 

18.2% 

11.4% 

96.7% 

58.0% 

20.0% 

12.7% 

— 35dB 

99.4% 

63.1% 

25.4% 

16.0% 

— 40dB 

99.5% 

68.0% 

30.6% 

18.5% 

Table  3:  Probability  of  false  target  detection  for  SNR 
interference.  Type  I  binary  mask  is  used. 


—20 dB  and  sinusoidal  FM 


SIR 

No  filtering 

STFT 

LPFT1 

LPFT2 

— 20dB 
— 25dB 
— 30dB 

28.5% 

17.0% 

2.1% 

2.0% 

75.4% 

18.0% 

3.7% 

2.7% 

96.3% 

26.4% 

5.2% 

2.7% 

— 35dB 

99.4% 

36.6% 

7.8% 

2.9% 

— 40dB 

99.4% 

43.5% 

10.4% 

5.1% 

Table  4:  Probability  of  false  target  detection  for  SNR 
interference.  Type  II  binary  mask  is  used. 


—20 dB  and  sinusoidal  FM 


SIR 

No  filtering 

STFT 

LPFT1 

— 20dB 
— 25dB 
— 30dB 

29.4% 

17.3% 

8.1% 

76.6% 

17.8% 

7.9% 

97.4% 

16.9% 

7.8% 

— 35dB 

99.6% 

19.3% 

8.7% 

— 40dB 

99.9% 

21.2% 

8.5% 

Table  5:  Probability  of  false  target  detection  for  SNR  =  —20 dB  and  LFM  interfer¬ 
ence.  Type  I  binary  mask  is  used. 


SIR 

No  filtering 

STFT 

LPFT1 

— 20dB 
— 25dB 
— 30dB 

29.5% 

2.9% 

1.5% 

76.2% 

6.4% 

1.8% 

96.4% 

8.1% 

1.9% 

— 35dB 

99.3% 

19.7% 

1.9% 

— 40dB 

99.7% 

45.3% 

2.7% 

Table  6:  Probability  of  false  target  detection  for  SNR  =  —20 dB  and  LFM  interfer¬ 
ence.  Type  II  binary  mask  is  used. 
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As  discussed  in  3.3  type  II  binary  mask  outperforms  type  I  binary  mask,  and  the 
LPFT  provides  better  results  than  the  STFT. 
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5  Conclusion 


The  problem  of  nonstationary  interference  suppression  in  noise  radar  systems  is  ad¬ 
dressed.  Towards  this  aim,  time- varying  filters  based  on  linear  TF  transforms,  namely 
STFT  and  LPFT,  are  developed.  The  filtering  is  performed  in  the  TF  domain  by 
using  a  binary  excision  mask,  which  removes  the  nonstationary  interference.  Two 
approaches  to  the  binary  mask  implementation  are  proposed.  Numerical  simulations 
show  that  the  TF  based  time-varying  filtering  significantly  improves  the  probability 
of  target  detection  in  severe  interference  environments.  The  best  results  are  obtained 
by  using  the  second  order  LPFT.  Results  clearly  demonstrate  that  the  time-frequency 
based  interference  suppression  method  can  significantly  improve  the  detection  per¬ 
formance  of  the  noise  radar  systems  in  extremely  strong  nonstationary  interference 
environments. 

The  proposed  interference  suppression  method  can  potentially  be  extended  to  other 
relevant  applications  such  as  over-the-horizon  radar  systems  (OTHR),  through- wall 
imaging  systems,  phased  array  radar  systems,  global  positioning  satellite  (GPS)  re¬ 
ceivers,  global  navigation  satellite  systems  (GNSSs),  etc.  Another  important  applic¬ 
ation  of  this  method  is  jammer  rejection  in  the  spread  spectrum  (SS)  communication 
systems.  Different  methods  have  been  proposed  for  rejection  or  mitigation  of  in¬ 
terferences  of  this  kind,  in  order  to  improve  interference  immunity  of  SS  systems 
and  provide  more  reliable  receiving  and  decoding  of  the  useful  signal.  The  proposed 
method  presents  another  viable  approach  to  jammer  mitigation  and  enhances  the 
performance  of  the  SS  receiver  in  such  severe  interfering  environment.  The  proposed 
method  may  also  be  extended  to  the  case  of  multiple  jammers.  This  method  can  gen¬ 
erally  be  applied  whenever  the  desired  signal  is  corrupted  by  broadband  interferences 
characterized  by  narrowband  instantaneous  bandwidths. 
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